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Abstract. The effects of gas accretion on spiral disk dynamics and stability are studied through N-body sim- 
ulations, including star formation and gas/stars mass exchange. The detailed processes of bar formation, bar 
destruction and bar re-formation are followed, while in the same time the disk to bulge ratio is varying. The 
accreted gas might be first prevented to flow inwards to the center by the bar gravity torques, which maintains it 
to the outer Lindblad resonance. While the first bar is weakening, the accreted gas replenishes the disk, increasing 
the disk-to-bulge ratio, and the disk self-gravity. A second bar is then unstable, with a higher pattern speed, 
due both to the increased mass, and shorter bar length. Three or four bar episodes have been followed over a 
Hubble time. Their strength is decreasing with time, while their pattern speed is increasing. Detailed balance 
of the angular momentum transfer and evolution can account for these processes. The gas recycled through star 
formation, and rejected through stellar mass loss plays also a role in the disk dynamics. Implications on the spiral 
galaxy dynamics and evolution along the Hubble sequence, and as a function of redshift are discussed. 
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1. Introduction 

Bars are an essential feature in galaxy evolution. About 
„ two-thirds of spiral galaxies are barred (de Vaucouleurs 

1963), one third being strongly barred (SB), the other 

third mildly barred (SAB). In the near-infrared, where 
' the central structure is unveiled from dust, the fraction is 

even higher (Eskridge et al. 2000). 

In the last decade, it has been realized that bars are 
not long lived features, in particular in gaseous spiral 
disks, and that galaxies are not frozen for a Hubble time 
in their morphological class: a bar can be destroyed by 
large radial gas inflow and mass accumulation in the cen- 
ter (Hasan & Norman 1990, Pfenniger & Norman 1990, 
Hasan et al. 1993, Friedli & Benz 1993). The mass concen- 
tration destroys the orbital structure that supported the 
bar, and this begins through the creation of two strong- 
inner Lindblad resonances (ILR). The process is initiated 
by the strong gravity torques exerted by the bar on the 
gaseous spiral arms, which decrease afterwards as the bar 
weakens. It can lead to the decoupling of a secondary bar, 
embedded into the primary, just inside its ILR (Friedli & 
Martinet 1993, Combes 1994), and this is likely to weaken 
the primary bar. 
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The bar phenomenon can then be a self-regulated pro- 
cess (e.g. Combes 2000). A tantalizing scenario is that a 
galaxy may have several bar episodes in its life, and that 
the fraction of barred galaxies in the de Vaucouleurs clas- 
sification only reflects the time percentage that a spiral 
galaxy spends as a barred object. How can a bar be re- 
vived after a destruction or weakening event? The disk 
self-gravity must be enhanced to counteract the stabiliz- 
ing influence of the bulge and central mass concentration, 
enforced by the first bar episode. This may be obtained 
through external gas accretion, that settles in the disk. 
Significant amounts of matter can be accreted in the life- 
time of the galaxy (e.g. Katz et al. 1996), and we consider 
realistic that the mass of a galaxy could double in a few 
Gyrs. The presence of gaseous warps in practically every 
spiral galaxy (Sancisi 1983, Briggs 1990) has been inter- 
preted as a sign that a spiral galaxy accretes mass and 
angular momentum: the amount accreted is such that it 
changes completely orientation in a typical time-scale of 
a few Gyrs (Jiang & Binney, 1999). 

During its life, a bar changes its pattern speed, and this 
could be a possible way to trace its age. In dissipationless 
simulations, the pattern speed decreases in a time-scale of 
a few Gyrs (Combes & Sanders 1981). This is also related 
to the length of the bar. The early bar instability involves 
central orbits, where the orbital precession rate is higher, 
then more and more orbits become trapped in the bar, 
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at larger radii, and the bar lengthens and slows down, 
The transient spiral waves formed in the stellar compo- 
nent transfer the angular momentum outwards. This pro- 
cess is also accelerated by escaping chaotic orbits, around 
corotation (Pfenniger & Friedli 1991). If there is a massive 
spheroidal dark matter halo, concentrated enough to per- 
turb the bar dynamics, the pattern speed of the bar could 
decrease even further, owing to the dynamical friction 
(Debattista & Sellwood 1998). According to the density 
of dark matter, the time scale could be estimated down 
to a few 10 8 yr. Since bars are observed most of the time 
to be fast rotators, this puts constraints on the density of 
dark matter in the central parts of galaxies. Fluctuations 
of the pattern speed can also be triggered by interacting 
companions (Gerin et al. 1990, Miwa & Noguchi 1998), 
but these remain transient. If gas accretion is able to reju- 
venate the bars, it is important to know the implications 
on the pattern speed, that could be a tracer of the event. 
This is studied in detail in the presented simulations. 

In Sect. 2 we briefly present the code and numerical 
methods used. In Sects. 3, 4 and 5 we present the results 
from our simulations. We give our conclusions in Sect. 6. 



In 2-D, the composition and the evolution of the bulge 
and the halo are not studied. They are only represented as 
analytical potentials. They are implemented in our sim- 
ulations as Plummer spheres. The associated potentials 
respectively read: 



(f>bulge(r) = G-M B 1 + 



and 



(t>haio{r) = G ■ M H ( 1 H — — 



(2) 



(3) 



The bulge and halo masses, Mb and Mjj, are varied ac- 
cording to table [l]. We here just note that Mb is gener- 
ally smaller than Mdisk, and that the halo mass included 
in the disk radius, M™^ r , which is only a fraction of 
Mb, is typically at most of the same order as Mdisk- The 
M B /M l d ™l lal and M H /M^l lal ratios are the two main 
parameters. On the other hand, the bulge and halo scale- 
lengths, as and a#, are fixed at as = 1-4 kpc and an = 25 
kpc (excepted in Sect. 3.3). 



2. Physical model for disk galaxies and numerical 
techniques 

To explore the influence of gas accretion on the disk dy- 
namics, there are many free parameters to vary; in a first 
step, to increase computing efficiency, we consider only 2- 
D models, with high spatial resolution and low CPU cost, 
with the constraint that the accretion has the same an- 
gular momentum direction as the initial galaxy disk. In a 
second step, with 3-D models, the angular momentum of 
the accreted matter is varied. 

We first describe the two-dimensional models and 
present the associated results in Sect. 3. We present 
the three-dimensional models, with a lower resolution, in 
Sect. 4. 



2.1. Disk, bulge and halo 

The galactic disk is first considered as an infinitely flat- 
tened two-dimensional component, embedded in a spher- 
ical halo, and a bulge partly stabilizes its inner regions. 
The initial disk is a Toomre disk, with a surface density 
of the form (Toomre, 1963): 



(i(r) = no 1 + 



'•disk 



(1) 



The scale-length ddisk has been varied between 4 and 6 
kpc. The disk is initially truncated at a radius Tdisk = 15 
kpc. This disk is made up of stellar and gaseous particles. 
Initially, the number of particles is 2 • 10 5 , the disk mass is 
M™* ia ' = 4 • 10 10 M , and its gas mass fraction is 30%. 
The disk mass Mdisk and the number of gaseous particles 
increase wit h time because of gas accretion, as described 
in Sect. 2.5. 



2.2. Gravitational interaction and gas dissipation 

The gravitational potential is computed via Fourier trans- 
forms. The particles are meshed on a useful Cartesian grid 
of 512x512 cells, and each particle mass is assigned to four 
grid points according to the " cloud- in-ccll" interpolation. 
Fourier transforms, carried out on a 1024 x 1024 grid in or- 
der to prevent Fourier periodic images do disturb the real 
disk, allow us to calculate the disk own potential, and de- 
rive forces. Then the analytical forces of the bulge and the 
halo are added, as they are described by (0) and (||). The 
potential differentiation and a bilinear interpolation give 
the gravitational force exerted on each particle. To sup- 
press two-body relaxation, inefficient in real galaxies due 
to the large number of particles, the gravitational forces 
are softened with a Plummer softening, of characteristic 
size of 200pc. 

To represent the dissipative character of the gas com- 
ponent, a frequently used method is the sticky-particles 
algorithm, which schematizes the cloud-cloud collisions, 
and the non-elasticity of the collisions, where the relative 
kinetic energy is lost to heating the gas, which efficiently 
re-radiates it away (e.g. Combes & Gerin 1985). The whole 
process tends to reduce the relative velocity of neighboring 
gas particles, cooling the gas through reducing its velocity 
dispersion. Here the process is represented by a dissipative 
force: 



/ = -k(v- < v >) 



(4) 



The local mean speed < v > is computed on a 1024 x 1024 
grid. The dissipation time-scale fc -1 is chosen between 5 
and 10 Gyr. 

The code finally computes the integration of equations 
of motion with a leap-frog algorithm (Hut et al. 1995). The 
time step for the gravity potential computation, and for 
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the integration of equations of motion, Atq 
time and is chosen to satisfy: 



At 6 < 



-/cell 



v m 



varies with 



(5) 



where L ce u is the cell size, and V max the maximum of the 
particles velocities. Computing V max allows us to make 
Atff, vary proportionally to L ce u/V max . The time step then 
varies between 0.4 and 1. Myr. Then the effective number 
of times steps to reach the Hubble time is about 20000. 

2.3. Initial velocities 

While creating the initial disk, with a spatial distribu- 
tion as described by (|l|), the code computes the mean 
circular velocity in each cell, the rotational frequency 
fl, the epicyclic frequency k, etc.. Random motions are 
given to all particles to prevent axisymmetric instabili- 
ties, their magnitude being at least equal to the critical 
speed (Toomre, 1964). The Toomre parameter Q is typ- 
ically Q = 1.5. The epicyclic approximation is used to 
relate the azimuthal and radial velocity dispersions: 



2n 



(6) 



The rotational velocity is reduced from the circular speed, 
according to the Jeans equations, applied to the infinitely 
thin disk. 



2.4. Star formation and stellar mass-loss 

The star formation scheme and the way the stars lose mass 
to the gas are inspired from those described by Jungwiert 
et al. (2001). The stellar mass-loss is not approximated as 
instantaneous, but decreases with time as a power-law. 

The star formation is assumed to be described by a 
local Schmidt law (Schmidt, 1959), i.e. the star formation 
rate is proportional to a power of the surface density of 
gas n£ as . We choose N = 1.4, according to the observa- 
tional determination of this parameter, on a global way 
and over large scales (see Kennicutt, 1998a, b). The gas 
density is computed on a 1024 x 1024 grid, then the star 
formation rate is calculated for each cell of this grid. In 
each cell, the mass of young stars is taken away from all 
gas particles in the cell, in proportion to their mass at this 
epoch. This new star mass is added to all stellar particles 
present in the cell, in proportion to their mass. If there 
is no stellar particle in the cell, we wait a few time steps, 
until such a particle comes into the cell. This delays the 
actual star formation, yet, for 95% of the disk surface (the 
whole disk excepted its outer edge), such a situation oc- 
curs less than once 20 time steps: so, the delay in the star 
formation is generally avoided. When some mass is redis- 
tributed from a gas particle to a stellar one, the associated 
linear momentum is redistributed, too. When a mass Am 
is transferred from particle 1 (gas) to particle 2 (stellar), 



with initial masses mi and rri2 and initial velocities v\ and 
V2, the masses and velocities of particles are changed to: 



mi 



Am 



m 2 = m2 + Am 



Do = 



777,2^2 + Amt>l 



(7) 
(8) 
(9) 



(10) 



This choice is not the only solution: the gas velocity 
is here conserved during star formation. One might also 
choose to add some velocity dispersion to new stars. 

The number of particles is then conserved, but the par- 
ticle masses vary. When a particle is accreted, it has the 
same mass as initially present gaseous and stellar parti- 
cles. The mass of stellar particles doubles in about 15 Gyr, 
while gaseous particles particles lose about half of their 
initial mass in 10 Gyr. 

Since each particle represents in fact already a cluster 
of actual components (star or gas cloud), the minimum 
mass of each particle being of the order of 10 6 M Q , the 
stellar mass-loss represents for each particle, the average 
mass-loss of an ensemble of stars, of various masses. As 
described in Jungwiert et al. (2001), we assume that a co- 
eval stellar population has a mass-loss rate M mass -i oss (t) 
equal to: 



M„ 



-ioss(^) — ^-initial 



Co 



t-t b + T 



(11) 



where Mi n m a i is the initial mass of the stellar population 
studied, and % its birth time. The two parameters are 
fixed to c = 5.35 • 10" 2 and T = 4.86 Myr. 

However, such a mass-loss description is suitable only 
for a population of stars born in the same stellar clus- 
ter. We may use it for the initial particles, but it is not 
convenient for particles containing stars born at different 
epochs. Yet we have to add the mass of young stars in 
particles containing older ones in the above star formation 
scheme. Such additions may change the mass-loss rate of 
particles. Thus, we will define the mass-loss law of any 
particle in a recursive way: 

— At t = 0, any particle has a mass-loss rate given by 
©■ 

— When new stars mass is added to a particle, we con- 
sider them as two merging particles: the old and young 
stellar particles. The first one is supposed to have a 
mass-loss rate given by: 



M old (t) = 



Cold 



t + T. 



old 



(12) 



The second one contains a coeval population, therefore 
it has a mass- loss rate described by (pJJ: 



t + T v , 



(13) 
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After the merging, the new particle has a mass-loss 
rate that can be written: 

Cn pin 



M new (i) 



t + T, n 



(14) 



We just have to determine the parameters c new and 
Tnew We impose the continuity of the total stellar 
mass-loss rate at the instant t of the fusion, which just 
reads: 



Cold 



t + T n 



t + T v 



t + T , 



■id 



(15) 



L young 

and the conservation of the total stellar mass-loss in- 
tegrated over 20 Gyr. These two equations determine 
the stellar mass-loss of a particle after the accretion of 
some young stars material. 

Once the stellar mass-loss is computed, we distribute 
this mass over the neighboring gas particles. We use a 
1024 x 1024 grid to compute the mass transfer in each 
cell. This mass is taken away from the star particles and 
transfered to the gas ones, in proportion to their mass, and 
the transfer of linear momentum that is associated with 
this mass redistribution is computed, just as was done for 
star formation. 

These processes are not taken into account at the same 
frequency as the particle advance, and gravity potential 
computations. The time steps for stellar mass-loss and star 
formation are: 



At ML = AtBF = 4Ai. 



2.5. Gas accretion 



(16) 



To represent gas accretion in the two-dimensional code, 
some assumptions are necessary. We assume that the 
mean speed of accreted particles is equal to that of the 
outer disk, considering that slower particles have been ac- 
creted earlier, and that the accreted gas random velocity 
corresponds to the same Toomre parameter as the disk 
(Q = 1.5). 

Accreted particles are placed just within the disk outer 
limit. None are put outside the disk. So, the disk radius 
is not directly changed in the process. The disk radius 
however increases indirectly, due to the disk response to 
accretion (see Sect. 3). The disk radius is then computed 
as a function of time, in order to continuously place new 
particles just below this radius, and not in the middle of 
the disk after some Gyr. 

The accretion rate, M accret iom is chosen to make the 
disk mass double in 7 to 10 Gyr (see for instance likely ac- 
cretion rates, derived by Jiang & Binney 1999). The gas is 
accreted from large-scale material. However, in some en- 
vironments, it is possible that the accreted mass depletes 
a significant fraction of the halo mass. Yet, we will note 
that the existence or not of a halo mass depletion does 
not influence the disk evolution. We then generally con- 
sider that the accreted mass is partly taken away from the 
halo (Mhaio — —Maccret/3), excepted in a few runs which 
are noticed in table [|. 



The consequences of all these assumptions are dis- 
cussed in Sect. [|. 

2.6. Three-dimensional code 

2.6.1. Disk, gravitation, and star formation 

The 3-D model is quite similar to the 2-D one, as far as the 
spatial distributions projected on the plane are concerned. 
In addition the mass distribution of the initial disk in the 
z direction is: 



2 I % 

p(r, z) — cr(r) sech — 
M 



(17) 



where the disk scale height is ho = 1.5 kpc, and a(r) is 
the disk surface density. 

The gravitational potential is still obtained by an FFT 
method. We now use a useful grid of 128 x 128 x 128. 
Proceeding Fourier transforms on a grid doubled in each 
dimension, as was done in the 2-D code, is not the best 
method as regards the computational efficiency. Instead 
we use the James (1977) method: appropriate masses are 
distributed on the grid boundaries to compensate for the 
Fourier periodic images. 

Star formation and stellar mass-loss are included in 
this code with the same schemes as in the 2-D one. 

2.6.2. Sticky-particles code 

The two-dimensional gas dissipation scheme cannot re- 
ally be used in three-dimensional simulations: computing 
the local mean speed in each cell implies that there are 
several gaseous particles in each cell. This would com- 
pel us either to use a few million particles, which means 
low computational efficiency, or to reduce the number of 
cells, which would make the spatial resolution reach a few 
kpc. Thus, we choose to use the sticky-particles algorithm 
(Schwarz, 1981). When two gaseous particles collide, their 
radial (along the line joining them) and tangential rela- 
tive velocities are reduced by factors (3 r and (it . We choose 
Pr = 0.65 and (3 t — 0.65 (see Jungwiert & Palous, 1996). 
Let us consider two particles at r\ and T2, with masses 
mi and 11J2, and velocities V\ and v 2 in the referential of 
their mass center . We note M = m\ + 7712, and: 



r 2 - n 



(18) 



Before the collision, the radial et tangential relative veloc- 
ities are respectively: 



v r = ((t>2 — Vi) ■ e r ) e r 

v t = (v 2 - vi) - ((t> 2 - vi) ■ e r ) e r 

After the collision, they are: 

v' t = PtVt 



(19) 
(20) 

(21) 
(22) 
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Then the new velocities of the particles in the referential 
of the mass center are: 

v'i = ^ (ftOi - v 2 ) + (P r - (i t ) {{ Vl - v 2 ) ■ e r ) e r )(23) 

If) -i 

v' 2 = -jj ((3 t (v2 - vi) + (J3r - A) ((«a - vi) ■ er) e r )(24) 

The energy lost after the collision is then: 

AE = (1 - (3 2 r ) v 2 r +(l- /3 t 2 ) v\ (25) 



We will discuss the influence of (3 r and (3 t on the results 
of our simulations in appendix 

2.6.3. gas accretion 

There are now many free parameters for gas accretion: 
the accreted gas angular momentum direction may be var- 
ied, and the accreted gas may be placed around any place 
of the disk, while in the 2-D code, it was always placed 
around the disk radius. We will precise which accretion 
parameters were used in each run when we study the 3-D 
results, in Sect. |[ 

3. Two-dimensional effects of gas accretion 

3.1. General observations 

3.1.1. Gas accretion and bar evolution 

Depending essentially on bulge to disk and halo to disk 
mass ratios, there are two kinds of bar evolution in 2-D 
simulations. We now describe them in turn. 

Bar maintenance Let us first examine the model of run 
Fl, the parameters of which are given in table [|. The disk 
evolution is plotted in Fig. [lj 

Without gas accretion, this disk develops a bar, with a 
short-lifetime. The bar is very strong around 1 and 2 Gyr. 
The central condensation rises very fast, and destroys the 
bar in less than 7 Gyr. The introduction of gas accretion 
makes the bar more long-lived. The disk has a bar wave 
during more than 15 Gyr. The bar destruction process still 
occurs, but is much slower. Therefore, one of the main 
effects of gas accretion is to maintain bars in such galactic 
disks. 

Moreover, the maintained bar can still be strong 
enough to correspond to a morphological SB-type, after 
more than 10 Gyr. A SB-type disk or strongly barred disk 
is a disk in which the bar is very long, with two arms con- 
nected to the bar. In a SAB-type disk, the bar is less long 
and has rather a lens aspect ; the spiral structure is rather 
disconnected from the bar. 

A strong bar is able to influence the spiral pattern 
outside of its corotation, so that the spiral and the bar 
have the same pattern speed. In the run Fl, at 15 Gyr, 
the spiral pattern appears always in phase with the bar, 
as in SB-type galaxies. Even if the bar is stronger during 
the first Gyrs, the disk may become more weakly barred, 



Run 


Bar lifetime 


Bl 


13.5 


B2 


15 


B3 


17 


B4 


12.5 


B5 


19 


B5h 


18.5 


B6 


13 


B7 


15.5 



Table 2. Bar lifetime for bar maintaining galaxies. 



as an SAB-type, and return to SB-type with an older bar. 
The frequencies of SB and SAB-type disks are studied in 



Sect. 3.3 



Runs B and Bl to B7 also show such a bar mainte- 
nance. The bar lifetime for each of them is described in 
table U 

Bar destruction and re-formation The second situation, ob- 
tained for instance with parameters of run F2 from table 
[j], corresponds to bar destruction and re-formation. The 
disk evolution of this run is shown in Fig. |^. 

Now, the isolated disk has a strong bar during less than 
2 Gyr. This bar destroys itself very fast. Then, a small 
lens-shaped bar remains in the inner disk, which testifies 
a previous strong bar existence. Even if the real bar is a 
short life-time one, the lens is observed over more than 10 
Gyr. 

Gas accretion allows the bar to reappear. After be- 
ing totally unbarred -even the lens is now destructed-, the 
disks undergoes a new bar episode. A cycle of bar de- 
struction and re-formation takes places in such disks. We 
generally obtain 2 or 3 cycles in 20 Gyr, sometimes 4. 

If the strength of the bar tends to decrease with time, 
because of the mass concentration, and the disk heat- 
ing (larger velocity dispersion), there still can exist quite 
strong bars (SB-types) after 10 Gyr. Due to gas accretion, 
and consequent disk cooling, the second or the third bars 
can be strong, almost as strong as the first bar. 

In runs A and Al to A9, we also follow such bar de- 
structions and re-formations. The number of bars followed 
over 20 Gyrs for these runs is given in table 

Dependence of bar evolution on parameters A galactic disk 
may either keep a single bar or get several successive bars. 
By simulating several disks, we study the dependence on 
physical parameters. We find that the two determining 
parameters are the bulge to disk and halo to disk masses 
ratios. Light bulge galaxies maintain their bar, while bars 
are destructed and regenerated in heavy bulge galaxies. 
The halo mass has also an effect on the disk evolution. 
The bar evolution of a galactic disk, depending on bulge 
and halo to disk masses ratios, is shown on Fig. |[ 

The bulge mass also influences the speed of the de- 
struction and re-formation cycle, in bar re-forming galax- 
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Run 
IXU.I1 


"disk 


O,halo 


1\ sinner 
M halo 


Mhulge 


M 

lvl accret 


M, , 

ivlhalo 


initial gas mass rate 


Fl 


4 kpc 


25 kpc 


2.5 10 10 Mq 


1.3 10 10 Mq 


0.6 10 10 Mq Gyr" 1 


accret / '2 


25 % 


0F1 


4 kpc 


25 kpc 


2.5 10 10 M 


1.3 10 10 Mq 





^[ accret /2 


25 % 


F2 


4 kpc 


25 kpc 


4.0 10 10 Mq 


4.0 10 10 Mq 


0.6 10 10 M Gyr" 1 


^[ accret /2 


25 % 


0F2 


4 kpc 


25 kpc 


4.0 10 10 Mq 


4.0 10 10 Mq 





^[ accret /2 


25 % 


A 


4 kpc 


25 kpc 


2.5 10 10 M 


3.0 10 10 Mq 


0.57 10 10 M Gyr" 1 


accret /3 


20 % 


OA 


4 kpc 


25 kpc 


2.5 10 10 M 


3.0 10 10 Mq 





^[ accret /3 


20 % 


B 


4 kpc 


25 kpc 


3.5 10 10 M 


1.4 10 10 Mq 


0.57 10 10 M Gyr" 1 


M-accret / 3 


20 % 


OB 


4 kpc 


25 kpc 


3.5 10 10 M 


1.4 10 10 Mq 





M^accret /3 


20 % 



Al 


4 


kpc 


25 


kpc 


2.0 


10 10 


Mq 


1.7 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


J^accret 


20 


% 


Alh 


4 


kpc 


25 


kpc 


2.0 


10 10 


Mq 


1.7 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 





20 


% 


A2 


4 


kpc 


25 


kpc 


2.0 


10 10 


Mq 


2.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / 3 


20 


% 


A3 


4 


kpc 


25 


kpc 


2.0 


10 10 


Mq 


3.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / 3 


20 


% 


A4 


4 


kpc 


25 


kpc 


3.0 


10 10 


Mq 


2.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret 


20 


% 


A4h 


4 


kpc 


25 


kpc 


3.0 


10 10 


Mq 


2.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 





20 


% 


A5 


4 


kpc 


25 


kpc 


3.0 


10 10 


Mq 


3.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr~ 


l 


accret / 3 


20 


% 


A6 


4 


kpc 


25 


kpc 


4.5 


10 10 


Mq 


2.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / 3 


20 


% 


A7 


4 


kpc 


25 


kpc 


4.5 


10 10 


Mq 


3.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / '3 


20 


% 


A8 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


2.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret 


20 


% 


A8h 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


2.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr~ 


l 





20 


% 


A9 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


3.5 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / 3 


20 


% 



Bl 


4 


kpc 


25 


kpc 


2.0 


10 10 


Mq 


1.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


M^accret / 3 


20 


% 


B2 


4 


kpc 


25 


kpc 


3.0 


10 10 


Mq 


1.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


M^accret / 3 


20 


% 


B3 


4 


kpc 


25 


kpc 


4.0 


10 10 


Mq 


1.7 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr~ 


l 


M^accret / 3 


20 


% 


B4 


4 


kpc 


25 


kpc 


4.5 


10 10 


Mq 


1.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


M^accret / 3 


20 


% 


B5 


4 


kpc 


25 


kpc 


4.5 


10 10 


Mq 


1.7 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


^^accret 


20 


% 


B5h 


4 


kpc 


25 


kpc 


4.5 


10 10 


Mq 


1.7 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 





20 


% 


B6 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


1.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / 3 


20 


% 


B7 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


1.7 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr" 


l 


accret / 3 


20 


% 



Rla 


4 


kpc 


25 


kpc 


2.5 


10 10 


Mq 


1.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


M^accret / 3 


20 % 


Rib 


4 


kpc 


25 


kpc 


4.0 


10 10 


Mq 


2.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


M^accret / 3 


20 % 


Rlc 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


3.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr" 


l 


M^accret / 3 


20 % 


R2a 


5 


kpc 


20 


kpc 


2.5 


10 10 


Mq 


1.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


M^accret / 3 


20 % 


R2b 


5 


kpc 


20 


kpc 


4.0 


10 10 


Mq 


2.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


M-accret / 3 


20 % 


R2c 


5 


kpc 


20 


kpc 


5.5 


10 10 


Mq 


3.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


— Maccret/3 


20 % 


R3a 


6 


kpc 


15 


kpc 


2.5 


10 10 


Mq 


1.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


— Maccret/3 


20 % 


R3b 


6 


kpc 


15 


kpc 


4.0 


10 10 


Mq 


2.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


— Maccret/'i 


20 % 


R3c 


6 


kpc 


15 


kpc 


5.5 


10 10 


Mq 


3.5 


10 10 


Mq 


0.6 


10 10 


Mq 


Gyr- 


l 


— Maccret/'i 


20 % 


Gl 


4 


kpc 


25 


kpc 


2.5 


10 10 


Mq 


2.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr- 


l 


M-accret / 3 


5 % 


G2 


4 


kpc 


25 


kpc 


2.5 


10 10 


Mq 


2.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr- 


l 


M-accret / 3 


15 % 


G3 


4 


kpc 


25 


kpc 


5.5 


10 10 


Mq 


2.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr- 


l 


M-accret / 3 


25 % 


G4 


4 


kpc 


25 


kpc 


2.5 


10 10 


Mq 


2.0 


10 10 


Mq 


0.5 


10 10 


Mq 


Gyr- 


l 


M-accret / 3 


40 % 



Table 1. Run parameters. The initial disk mass is always M^"^ a/ = 4 • 10 10 Mq. The runs Fl and F2 are introduced 



in Sect. 3.1. The runs A and B are studied in Sect. 3.2. Runs Al to A9 and Bl to B7 are used to determine the 



dependence of bar evolution on parameters ; in terms of bar evolution, these runs are rather similar either to run A 



or to run B. Runs Rla to R3c allow us to study the link with the de Vaucouleurs classification (see Sect. 3.3) 



Fig. 1. Runs Fl and 0F1: disk evolution with accretion and bar maintenance (Fl, top), and without accretion (0F1, 
bottom). 
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Fig. 2. Runs F2 and 0F2: disk evolution with accretion and bar destruction and re- formation (F2, top), and without 
accretion (0F2, bottom). Here, the bar is totally destroyed before its re-formation. It can also be much weakened 
without disappearing (see Fig. |i| for such a situation in model A). Two bar episodes can be seen here in 20 Gyrs ; 
three or four bar episodes have been followed as well in similar runs (A1-A9). 



Run Number of bars in 20 Gyrs 



Al 


4 


Alh 


4 


A2 


2 


A3 


2 


A4 


3 


A4h 


3 


A5 


2 


A6 


3 


A7 


2 


A8 


3 


A8h 


3 


A9 


2 



Table 3. Number of successive bars followed over 20 Gyrs 
in bar re-forming galaxies. 



A A 


♦ 




♦ 




[13.5] [15.5] I 


(3) 




(2) 




A A 


♦ 




♦ 




[12.5] [19] / 


(3) 




(2) 


■ F2. 










(2) 


TB / 










A [1 % / 


♦ 




♦ 




[15] [17] / 


(3) 




(2) 




TF1 / 




■ A 






a m A 


♦ 


(3) 


♦ 




[13.5L- (A) 
i " i 


(2) 


i 


(2) 


i 



1.5 



0.5 



0.25 



0.5 



0.75 
M bulge /M disk 



A Runs B1-B7: Bar maintenance, 
[bar lifetime in Gyr] 

♦ Runs A1-A9: Bar destruction and re-formation, 
(number of bars over 20 Gyr) 

Fig. 3. Disk evolution depending on bulge to disk and 
halo to disk masses ratio. Mdisk is the initial disk mass, 
and M l h n Z r is one tenth M halo . 



ies. It influences the bar lifetime in bar maintaining galax- 
ies. Both influences are explained on Fig. ||. 

The initial gas to stellar disk mass ratio has also been 
varied (runs G1-G4). The higher it is, the stronger the 
density waves are during the 1 or 2 first Gyr. Yet, this 
ratio does not really influence the long-term evolution of 
the disk: After a few Gyr, the gas to stellar mass ratio 
depends more on the accretion rate rather than on the 
initial gas to stellar mass ratio. 



3.1.2. Other consequences of accretion 

The main other consequences of accretion are: 
— the spiral structure is rejuvenated by accretion (see 



also Sect. 3.3). Arms have roughly the same morphol- 
ogy around 15 Gyr as around 3 Gyr, whereas in non 
accreting disks, they become weaker, since the gas den- 
sity then decreases in the disk while the gas accumu- 
lates in an inner ring and/or an outer ring, 
the disk radius increases; accreted gas, even if always 
placed inside the disk border, produces larger disks, 
there is sometimes a gaseous outer ring, which is 
stronger than without accretion. This, as shown be- 
low, is due to the accreted gas repelled outside the 
corotation because of bar torques, 
the bar pattern speed raises; this will be shown and 
interpreted in Sect. 3.2. 



The three last phenomena are now studied and explained. 

3.2. Detailed study of two galactic disks 

Let us now present two runs that have been studied in 
detail, runs A and B in table [lj 

Model A corresponds to a disk which undergoes three 
successive bars, whereas model B corresponds to a single 
bar that persists over a Hubble time. The bar strength 
evolutions for both models are shown in Fig. |J. The bar 
strength is defined as the ratio of the tangential force of 
the to = 2 Fourier component, normalized to the radial 
force. If the potential is decomposed as 

*(r, 6) = * (r) + cos ( md ~ <t>™) (26) 

m 

the bar strength is the maximum over radius of 



S 2 = 2$ 2 /r|F r 



(27) 



where F r = —d^o/di. 

The cause of these two different scenarii in models A 
and B, the evolution of the bar pattern speed, and the 
increase in the disk radius, are now developed and studied. 

3.2.1. Evolution of the accreted gas distribution 

The accreted gas is first located in the outer disk. Yet, 
non axisymmetric gravity fields and dissipation make its 
distribution evolve. The accreted gas distribution and its 
evolution are shown in Fig. |5|. Only particles that are ac- 
creted during the first 1.5 Gyr are taken into account in 
this figure, in order to have a reasonable time axis, as the 
accretion period of such particles is rather short in com- 
parison with the galaxy evolution time. If we consider the 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 

time (Gyr) time (Gyr) 



Fig. 4. Bar strength evolution in models A (left) and B (right). 

Fig. 5. Evolution of the radial distribution of accreted during the first 1.5 Gyr gas. Model A is on the left, model B 
is on the right. 




Fig. 6. Second harmonic of the potential, as a function of radius and pulsation, for model B. The main maximum 
(at small radii and high frequencies) corresponds to the bar, the secondary one (at large radii and low frequencies) to 
spiral waves. Similar analysis are done for model A. 



particles that are accreted between any epoch t and t+1.5 
Gyr, the results are very similar to those plotted in Fig. ||. 

For both models A and B, we see that the accreted gas 
density is never high in the region between 4 and 7 kpc, 
which corresponds to the bar corotation in both models. 
Tangential forces exerted by the bar prevent gas from stay- 
ing in this part of the disk. Then, accreted particles may 
either stay in the outer disk, or fall towards the disk center, 
if some dissipation make them go there. Actually, we see 
on Fig. [5] that accreted gas accumulates around the galac- 
tic center after some Gyr. To go there, it must overcome 
bar positive torques. The only mechanism responsible is 
the gas dissipation, which becomes higher and higher as 
gas density outside the corotation increases. By spreading 
out in the disk, some gas reaches the corotation. Then, 
bar negative torques drive it quickly into the disk center. 



Yet, there are some differences between models A and 
B. The gas repelled outside the corotation seems to be 
more abundant for model A. The particles get over the 
corotation later, after 8 Gyr instead of 6 or 7 Gyr for 
model B. This interval is rather small, but it is signifi- 
cant that for model B, there is still a strong bar when the 
gas falls inside the corotation, whereas in model A, the 
corotation crossing occurs only when the bar strength is 
minimum. Moreover, the gas infall inside the corotation 
is rather regular and takes some Gyrs for model B ; for 
model A this infall is faster and happens at once, between 
the two first bar episodes. 

This difference concerning the accreted gas repartition 
allows us to account for the two encountered situations. 
Moreover, it will explain that, in the next subsections, 
different phenomena occur in a continuous way for model 
B, and only between bar episodes for model A. 
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time (Gyr) time (Gyr) 



Fig. 7. Bar pattern speed as a function of time, for model 
A (full line) and model B (dashed line). 

3.2.2. Explanation of bar destruction and re-formation 
or maintenance 

The common point between models A and B is that gas ac- 
cretion favours bars. The reason for that is the increase in 
the disk mass. In non accreting disks, the central conden- 
sation mass raises, which stabilizes the disk against barred 
waves. Accretion makes the disk mass increase and become 
more self-gravitating, while making the central (inner disk 
and bulge) become relatively lower, which may either pre- 
vent the bar from being destroyed or re-form a bar. 

The difference takes place in the accreted gas repar- 
tition. For model B, this gas falls inside the corotation 
continuously. So, the region of the possible bar (which is 
r < 4 to 6 kpc) increases in weight regularly. The bar is 
then favored regularly, and is maintained a long time. 

For model A, the effects favourable to bars occur only 
between the bar episodes. During such a bar episode, the 
accreted gas stays in the outer disk. Even if the total disk 
mass rises, only the outer regions increase in mass: this 
cannot maintain the bar nor create a new one. Between 
two bar episodes, the region of the possible bar gains mass. 
The possible barred region to central condensation mass 
ratio becomes lower, and a new bar wave appears. 

Let us consider now the bar pattern speed evolution, 
and how it is related to the accreted gas distribution. 

3.2.3. Influence of accretion on bar pattern speed 

The bar pattern speed, Of, , has been measured via tem- 
poral Fourier transform of the m = 2 harmonic of the 
potential. Fig. ^| shows the intensity of this harmonic as 
a function of radius and frequency ft, around 3.5, 9.5 and 
15.5 Gyr. Using such Fourier analysis, we are able to sum- 
marize ftb as a function of time, as shown in Fig. 0. 

This reveals a significant increase in the bar pattern 
speed, which may double over a Hubble time. The differ- 
ence between the two A and B cases is that this growth 
is either smooth or irregular. This suggests that the bar 
acceleration is closely linked to gas accretion, since it has 



Fig. 8. Evolution of the disk inner mass, Mdi S fe(2kpc), 
normalized to its initial value, for model A (full line) and 
model B (dashed line). Bulge and halo inner masses are 
to be added to it to give Mi nner . 

the same evolution as the accreted gas distribution in both 
models A and B. 

The implications of such a pattern speed increase are 
developed below. Let us first analyze the reasons for this 
phenomenon. 

3.2.4. Physical cause of the bar acceleration 

Central condensation growth A possible cause of the 
growth in the bar pattern speed is the inner mass ac- 
cumulation. For every particles, this would make ft and 
k rise. If we write M(r) the total (disk, bulge and halo) 
mass that is inside the r radius, all frequencies for par- 
ticles around r arc to increase proportionally to y/ ' M{r). 
We will study the evolution of Mi nner = M(2 kpc), for 2 
kpc is generally about the mean radius of the bar compo- 
nents (whereas the maximum of the potential harmonics 
are around the bar extremity). 

We obtain the expected increase in Mi nner (see Fig. |J). 
The aspect of this growth, regular or not, proves that there 
is a link between it and bars succession. For model A, each 
bar episode makes M inner raise: between two bars, much 
accreted mass comes inside the corotation, then the next 
bar exerts tangential forces that make this mass fall into 
the central condensation. For model B, Mi„ ner is rather 
constant: accreted gas crosses the corotation in a perma- 
nent way, and supplies the central condensation during 
the whole galactic evolution. Let us note that Mi nner is 
considered as constant (for model B), while it is actually 
smaller during the 6 or 7 first Gyrs. Yet, this is explained 
by the absence of particles accreted before t = 0: most of 
the actually accreted particles only come into the inner 
disk after 6—7 Gyr, and only few of them go into the 
central condensation earlier. So, this change in Mi nner is 
linked to the beginning of accretion at t = 0, it has no 
physical meaning, and what is significant is that Mi nner 
is constant for t > 7 Gyr. 
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growth factor between the first bar and... 


\/ ^"firmer 




...the second bar 


1.06 


1.11 


...the third bar 


1.20 


1.37 



Table 4. Inner mass and bar pattern speed growth rates 
(model A). 
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16 



Fig. 9. Inner (r < 2 kpc) stellar and gaseous particles 
angular momentum per unit mass, normalized to its value 
at t = 0, for model A (full line) and model B (dashed line). 



Accretion is responsible for an increase in Mi nner . But 
this increase is not sufficient to explain the whole bar ac- 
celeration. For instance, the growth rates of V Mi nner and 
of fij for model A, between several instants, are given in 
table |^. The increase rate of \/Mi nner is not sufficient, 
and is only partly responsible for the increase in fit,: the 
growth rate Of, is twice bigger. The bar pattern speed is 
also influenced by another phenomenon. 

Accretion and angular momentum supply In non accreting 
galactic disks, non axisymmctric fields provoke angular 
momentum transfers from bar to gas outside the corota- 
tion, and from gas inside the corotation to bar. In model 
A, the gas density is high outside the corotation when 
there is a strong bar, then the bar is expected to give some 
momentum to outer particles. When the bar strength is 
lower, much gas arrives in the inner disk, and may bring 
there its angular momentum. This is observed on Fig. ^|. 
For model B, the accreted gas density evolves more reg- 
ularly, and the angular momentum per unit mass raises 
steadily. 

Gas accretion has therefore the main implication to 
raise the inner angular momentum per unit mass. If £1 is 
higher for a particle, its orbit's precession rate will also be 
faster, and it is likely that the bar instability it will par- 
ticipate in, will have a higher pattern speed. The angular 
momentum supply makes the bar pattern speed raise. 

Moreover, gaining angular momentum should make 
the bar become more circular, because faster particles have 
less eccentric orbits. This is observed on Fig. [To] for model 
A, and this phenomenon is similar for model B. 



5 20 




10 

time (Gyr) 

Fig. 11. Disk radius as a function of time, for models A 
(full line) and B (dashed line). 

The main implications of gas accretion are therefore an 
increase of the central mass and of the angular momentum 
supply, and a higher bar pattern speed. Moreover, the bar 
becomes rounder as it gains angular momentum. 

3.2.5. Evolution of the disk radius 

The last point to discuss about models A and B is the disk 
radius augmentation. This is plotted on Fig. |ll]. Here, the 
disk radius is defined as the radius containing 90 perccnts 
of the disk mass. 

The stronger the bar, the faster the radius grows. For 
both models, the OLR is outside the disk. Yet, bar torques 
may not explain totally the radius growth, for at least two 
reasons: 

— even if they could enlarge outer particles orbits, they 
would have the same consequences without any gas 
accretion, which is not observed. 

— when models are studied in which the OLR is inside the 
disk, the bar torques could only make the disk radius 
decrease. Yet, we still observe a radius increase. 

The likely explanation is rather that the stronger the bar, 
the more accreted gas is accumulated in the outer disk: 
this gas spreads through the disk because of its dissipation; 
a part goes inside the corotation, and a part makes the 
disk expand. 

3.2.6. Bar pattern speed and corotation radius 

We have shown that gas accretion makes the bar pattern 
speed change. Yet, this does not imply that any change 
could be easily observed in real disks: the bar pattern 
speed cannot be directly observed, and the corotation ra- 
dius is measured instead. So, we may wonder how the 
corotation radius evolves, and compare it to the disk ra- 
dius. For instance for model B, the corotation radius and 
the corotation to disk radius ratio are given on Fig. |lj. 

The corotation radius decreases with time (see 
Fig. |l2|) . This can be interpreted in terms of bar pattern 
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Fig. 10. Inner disk gaseous and stellar density for model A, seen at 3, 10 and 15 Gyr. The unit of length is 300 pc. 

0.4 




cc 

o 
o 
CC 



10 12 
time (Gyr) 




10 12 
time (Gyr) 



Fig. 12. Corotation radius (left) and corotation to disk radius ratio (right) for model B. 



speed: we have noticed that the pattern speed fib raises 
with time faster than the square root of the inner mass, 
thus faster than the particles angular velocity £l(r). Then, 
the corotation radius R CO r, defined by = £l(R cor ), is to 
decrease with time. Moreover, we know that the disk ra- 
dius becomes larger, so that the corotation to disk radius 
ratio, which is the more relevant parameter, decreases in 
a significant way as displayed on Fig. [l2j 

The bar pattern speed evolution is related to a change 
in the corotation radius, and makes the corotation to disk 
radius ratio decrease with time. Moreover, even if the bar 
pattern speed does not change a lot during the first Gyrs 
(see Fig. ^), the change in the corotation radius is larger 
during the first Gyrs (see Fig. |l2|); then the position of the 
corotation is significantly modified by gas accretion even 
in rather young disks. 

3.3. Gas accretion and spiral morphology 

We have already briefly noticed that gas accretion rejuve- 
nates spiral patterns. In an accreting disk, the main struc- 
ture of the outer disk consists in spiral arms, whereas in 
non accreting disks, such a structure is rather evanescent: 
it is replaced by circular structures in about 8 Gyr. 

Actually, the gas mass fraction is obviously higher; so 
the stellar and gaseous medium are globally cooler: this 
prevents the stability parameter Q from increasing too 
much, so that the spiral structure persists in the disk. 
When there is no accretion, the first density waves heat 
the medium, which does not cool back fast enough because 
of lack of gas. 

Moreover, the accreted gas spreading prevents the disk 
from getting almost empty of gas in specific regions (even 
if does not prevent rings from appearing, as will be shown 
later). The initial gaseous density radial distribution is 
better preserved, except a change in its scale length. This 
let the gaseous medium stay more favourable to spiral 
waves propagation. 



Rejuvenating the spiral structure does not prevent res- 
onant ring formation. Such rings may even be stronger. 
External rings are intensified by accreted gas accumula- 
tion in the outer disk. Inner rings, connected to ILR, do 
not seem to be affected by accretion. The most favored 
rings are those that surround the bar: they are related 
to the 4/1 resonance. They are weaker than other rings 
in non accreting disks, due to gas depletion near coro- 
tation. Accretion lets them become stronger. The subse- 
quent disks are then typically RSX, SX(r), or even RSX(r) 
(X meaning A, B or AB) in the de Vaucouleurs classifica- 
tion. In this classification, the meaning of SA , SAB and 
SB has already been explained (see Sect. 3.1). In a S A 
disk, the residual bar strength, as defined in Sect. 3.2 is 



less than 0.2, In SAB disks, the bar strength is contained 
between 0.2 and 0.5 ; and in SB disks, it is more than 0.2 
and can be much larger. We notice that for moderate bar 
strengths (0.3-0.4), the disk can be either SAB or SB. A 
RSX-type disk shows an outer ring. The relations between 
bar strength and disk morphology have been studied by 
Block ct al. (2001). A SX(r)-type disk shows a ring sur- 
rounding the bar: it is not an inner ring connected with 
ILR, but a larger ring. This SX(r)-type is much favored 
by accretion. 

Let us try to interpret the models in terms of this clas- 
sification, exactly with the distinction between SB and 
SAB galaxies. In SB disks, spiral arms generally appear 
in phase with the bar and might have the same pattern 
speed; in the weaker SAB ones, spiral patterns are some- 
times separated from the bar, and are likely to have a 
lower pattern speed. 



The two disks studied in Sect. 3.2 can be classified SA 
or SAB, according to this definition, since Fig. 



reveals 

an angular speed lower for arms than for bars. However, 
the disk sometimes presents an SB aspect, for both models 
A and B. But the frequency of strong bars (SB) does not 
exceed 30%, and the rate of SB disks in the whole galactic 
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Run 


SB to [SAB+SB] ratio 


SB global rate 


Run 


M bulga /M%T l 


Rla-b-c 


20 to 30 % 


15 to 20 % 


1 


0.8 


R2a-b-c 


40 to 45 % 


25 to 30 % 


2 


0.6 


R3a-b-c 


75 to 85 % 


50 to 55 % 


3 


0.3 



Table 5. Rate of SB-type disks for different parameters. 



evolution is less than 20%, whereas observed galaxies are 
divided equally between the three categories, thus half of 
barred galaxies are of SB-type. 

In fact, the frequency of SB disks in simulations de- 
pends on some input parameters. The halo to disk and 
bulge to disk ratios do not intervene there, neither does 
O'buige- The important parameters are cidisk and 

Clhalo- It 

is to be noted however that when ahaio is varied, Mhaio is 
varied also in order to keep M^™" er constant, which is the 
relevant mass. Moreover, increasing ahaio is equivalent to 
decrease cidisk- The runs Rla to R3c are then used here, 
and the frequency of SB disks as a function of both scale- 
lengths is given in table ||. Approximatively half of barred 
disks are of SB-type, and a third of disks are SA, a third 
SAB, and a third SB, which is in good agreement with 
observations. 

Does accretion have some influence on the SB to SAB 
ratio (obviously, it reduces the SA rate, and raises the 
SAB and SB global rate, for it favours bars)? If we stop 
accretion at any epoch of simulation, the SB frequency 
decreases; furthermore, non accreting disks have lower SB 
frequencies than accreting ones. So, gas accretion favours 
a high SB to SAB ratio. Accretion brings mass in the outer 
disk, which makes adisk rise, whereas in non accreting 
disks, the bar makes a heavy central condensation, which 
on the contrary reduces the disk scale-length. Finally, gas 
accretion favours not only bars, but also strong SB bars. 

4. Three-dimensional effects of gas accretion 

We now present the results of 3-D simulations concerning 
bar evolution. New free parameters arise, in particular the 
angle between the angular momentum of the accreted gas 
and that of the disk. Also, the gas can be accreted at any 
radius. For the sake of comparison with the previous 2-D- 
results, let us begin by in-plane accretion. 

4.1. In-plane accretion 

In these first runs, as in the 2-D-case, the gas is assumed 
to be accreted in the plane, at all azimuths in the disk, 
with the tangential velocity corresponding to the rotation 
curve. This situation corresponds to a long-term, quasi- 
stationary accretion over a Hubble time, as if the optical 
disk was surrounded with slowly infalling gas, that has 
already settled into a larger-scale disk. Parameters of the 
runs are given in table [| 

For runs 1 and 2, we observe bar destruction and re- 
formation. These runs show respectively 2 and 3 successive 



Table 6. Parameters of the 3-D runs with in-plane accre- 
tion. The value of Mf™f° r /M d ^ al is fixed at 0.75. The 
scale-length for the disk is 5kpc, for the bulge 1 kpc, and 
for the halo 20 kpc. The accretion rate makes the disk 
mass double in 7 Gyrs, and the initial gas fraction is 20%. 



bars over 15 Gyr. Run 3 presents a long-lifetime bar, that 
is still present in the disk after 13 Gyr. The bar accelerates 
such that its pattern speed doubles in about a Hubble 
time. The accreted gas distribution in the disk, related to 
bar destruction and re- formation in runs 1 and 2, is similar 
to that obtained in 2-D models. The vertical distribution 
of particles is not modified by gas accretion in these runs. 
Therefore, the results and interpretations reported for 2-D 
simulations are strengthened by these new runs. 

4.2. Out-of-plane accretion 

When gas is not accreted in the galaxy plane, we define 
Race the ratio between the accreted gas mean radius and 
the disk radius (see Fig. [b|). The angle between the ac- 
creted gas angular momentum and the initial disk angular 
momentum is V'acc- The accreted gas mean initial veloc- 
ity has again no radial component. It seems reasonable to 
assume that R aC c and ip acc are nearly constant over a few 
Gyrs, but not over a Hubble time. Therefore, the accre- 
tion phase lasts only 5 Gyrs. It is started well after the 
first relaxation phase, when a first bar develops in the ini- 
tial disk ; a few Gyrs later, the bar weakens, and without 
gas accretion, the galaxy would not be strongly barred 
any more. We hereafter note t = the beginning of the 
accretion event, and not the beginning of the whole disk 
evolution. The accretion does not take place from all az- 
imuths around the disk (at radii around R aC c), but only in 
a 3-kpc-sized area, which seems more realistic for out-of- 
plane accretion. The mass accreted over 5 Gyrs is half of 
the initial disk mass. The disk evolution is followed during 
the accretion episode and later. Two important parame- 
ters, V'acc and Race, are varied. Parameters of the run are 
given in table 0. 

4.2.1. Large-radius accretion 

In runs ACC8-0 to ACC8-87, the gas is accreted in the 
outer disk with several different values of R acc . The gas 
accretion always re-inforces the bar significantly. Warps 
are also formed in the outer disks, while the initial disk 
axis rotates. 

General disk evolution When gas is accreted with an angu- 
lar momentum inclined to that of the disk, it is submitted 



Bournaud & Combes: Gas accretion on spiral galaxies: bar formation and renewal 



13 



Accreted gas 




Fig. 13. Model for gas accretion and parameters R aC c and 



Run 


Race 




M bulge /Msi:r i 


REF 




no accretion 


0.8 


ACC8-0 


0.8 





0.8 


ACC8-8 


0.8 


8 


0.8 


ACC8-20 


0.8 


20 


0.8 


ACC8-32 


0.8 


32 


0.8 


ACC8-45 


0.8 


45 


0.8 


ACC8-65 


0.8 


(>") 


0.8 


ACC8-87 


0.8 


87 


0.8 


ACC8-20-2 


0.8 


20 


0.6 


ACC8-20-3 


0.8 


20 


1.0 


ACC9-20 


0.9 


20 


0.8 


ACC7-20 


0.7 


20 


0.3 


ACC6-20 


0.6 


20 


0.3 


ACC3-25 


0.3 


20 


0.8 



Table 7. Parameters varied for out-of-plane accretion. 
R acc is the ratio between the accreted gas radius and the 
disk radius, and ip&cc is the angle between the angular mo- 
menta of the accreted gas and the initial disk, as described 
on Fig. ||. We choose M^ r /M^ al = 0.75. 



Fig. 14. Edge-on views of the gaseous component at t=6 
Gyr (1 Gyr after the end of accretion) for run ACC8-45. 
The view axis is perpendicular to the major axis of the 
bar (top), then parallel to it (bottom). For this run, the 
change in disk orientation is almost over at t=6 Gyr. 

to restoring forces from the disk, that make the orbits pre- 
cess. The precession is differential, and the various gaseous 
orbits, at different radii will soon intersect, and dissipa- 
tion will occur and make newly accreted gas align to the 
disk plane. Therefore all the matter will gather in a disk, 
which axis will slowly change orientation in a few Gyrs 
(see Jiang & Binney, 1999). This change in the disk orien- 
tation is observed in the present simulations (see Fig. |l4| ). 
There is also a flaring in the re-oriented gaseous disk, while 
this flaring disappears during the change of orientation. A 
polar ring also appears in run ACC8-87. 

Bar strength and out-of-plane accretion For runs ACC8-8 
to ACC8-87, a bar appears in the disk. It is very long 



Fig. 15. Stellar and gaseous particles viewed perpendicu- 
lar to the re-oriented disk at t= 5 Gyr, i.e. the end of the 
accretion event, for run ACC8-45. This disk shows two 
spiral arms in phase with a very strong bar. 

around t=3-5 Gyr, as shown on Fig. [r] for run ACC8- 
45 at t=5 Gyr. For this runs, the maximal bar strength 
is 0.75, and exceeds 0.85 for runs ACC8-0, 8 and 20. In 
the non- accreting disk REF, the primordial bar is weak- 
ened at t — and the bar strength is never more than 
0.3. Bars re-formed by out-of-plane gas accretion are much 
stronger than the bar re-formed in the in-plane accretion 
runs, where the bar strength never exceeded 0.5-0.6 in 
either 2-D or 3-D runs, whatever the different parame- 
ters were. The lower strength of previous bars are obvious 
when comparing Figs. |l| and § 

Thus the present out-of-plane gas accretion model ap- 
pears more favourable to barred waves than the previous 
in-plane accretion models. To interpret this, let us consider 
the run ACC8-0: this disk also develops a very strong bar, 
of maximal strength 0.9, while the angle of accretion is 
zero. In fact, the common new feature introduced in the 
3-D runs is to limit the accretion to a 3-kpc-sized area 
of the disk. This increases the non-axisymmetric compo- 
nent of the gravitational potential and favors a stronger 
bar. Moreover, when ip acc > 0, the accretion paths inter- 
sect the disk plane around two points: this strengthens 
m = 2 harmonic of the potential, which may explain the 
enhancement of two spiral arms in the disks. 

In summary two phenomena may explain why gas ac- 
cretion favor bars in disk galaxies: 

— the supply of mass to the disk, which reduces the rel- 
ative mass of the central stabilizing condensation. 

— the asymmetric distribution of the accreted gas. 

It is difficult to disentangle the dependence of each of them 
on the accretion parameter V'acc, since both phenomena 
occur simultaneously. Yet, the main trends are: 

— with larger i/'acc, the gas accretion onto the disk is 
longer and slower, which does not favor bar formation. 

— with larger ^ acc , the distribution of the accreted gas 
is more asymmetric, and the m = 2 harmonic of the 
resulting potential is higher. 

The two mechanisms to form, re-form, or strengthen a 
bar, have thus opposite dependence on V'acc- The evolu- 
tion of the actual bar strength on ip &cc will allow us to 
determine which is prevailing. The bar maximal strength 
as a function of ip acc for runs ACC8-0 to 45 is plotted on 
Fig. [l§ 

We find that the maximal bar strength decreases when 
ip aC c increases. It appears that the influence of a strong 
and rapid mass supply dominates the process of bar re- 
formation. 

Lifetime of the bar and gas accretion Let us now examine 
the dependence of the bar lifetime on ip aC c- Considering 
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Fig. 16. Maximal bar strength over the whole bar evolu- 
tion, as a function of ipa.cc, with R acc =0.8. 
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Fig. 18. Bar pattern speed as a function of time for several 
runs. 
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Fig. 17. Lifetime of the bar as a function of the accre- 
tion parameter ip acc . The bar episode ends when the bar 
strength is less than 0.2. The mean accretion radius is 

Racc=0.8 

that the bar is destroyed when its strength drops below 
0.2, the lifetime of the bar in each run is shown on Fig. [H]. 
The first point is that all these bars are rather short-lived 
features. This is not surprising, for at least two reasons: 
first the gas accretion is supposed to be limited in time, 
and occur during only a few Gyrs; second the formed 
bars are very strong, so the induced mass transfers de- 
stroy them in only a few Gyrs. The second point is the 
dependence of the bar lifetime with the accretion angle. 
The bar lifetime is longer for smaller V'acc, and therefore 
the strongest bars have here the longest lifetimes. It is 
well-known that in non-accreting disks, the strongest bars 
are the most evanescent ; thus the bar lifetime is strongly 
influenced by gas accretion. A low ip acc favors the mass 
supply to bars, which lengthens their lifetime. Accretion 
not only creates a bar, but also maintains it. 

Evolution of the bar pattern speed The evolution of the 
bar pattern speed with time is shown for different runs on 
Fig. [l8|. Two steps are clearly seen: when the bar is rather 
young (i.e. before 2-3 Gyrs), the pattern speed rises. Then, 



after a peaked maximum, it decreases while the bar is de- 
stroyed. Without accretion, the bar rotation would always 
slow down (Combes & Sanders, 1981) ; the acceleration 
observed here at the beginning is related to gas accretion. 
This succession of acceleration and deceleration of the bar 
rotation is similar to that observed for 2-D models with 
bar destruction and re-formation (see Fig. Q). 

In summary, out-of-plane gas accretion, like in-plane 
accretion, forces the bar pattern speed to increase, while 
bars naturally tend to slow down. 

4.2.2. Small-radius accretion and dependence on 
accretion radius 

When gas is accreted from outside the plane, the choice of 
accretion radius R acc is free. In the simulation ACC3-25, 
the accretion radius is chosen very small. The disk radius 
decreases in a few Gyrs, and no bar is formed in the disk. 
However such a situation is rather unusual, since gas ar- 
riving at larger radii will be accelerated in the infall, and 
the accretion will take a long time. A rapid deceleration 
requires violent conditions and strong dissipation, that oc- 
cur in galaxy collisions and mergers. A very small radius 
of accretion with high mass accretion rate is not likely to 
resemble a quasi-steady accretion. 

Let us now study the influence of the accretion radius 
when R acc is of the same order as r^isk '■ the runs involved 
are ACC6-20, 7-20, 8-20 and 9-20. There is no significant 
dependence of the bar evolution on R acc . In fact, in the 
above runs, a weakened lens-shaped bar is present in the 
disk when we start accreting gas: its corotation radius is 7 
kpc. Then, for runs ACC6-20, 7-20, 8-20 and 9-20, the gas 
is accreted outside the corotation (the accretion radius is 
more than 7.5 kpc in those runs). As we know from 2-D 
studies that the important factor is the relative position of 
the accreted gas with respect to corotation radius, we can 
understand that the accretion radius does not influence 
the bar evolution, provided that this accretion radius is 
larger than the corotation one. 
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4.3. General situations 

There are many free parameters concerning the way gas 
is accreted on a galactic disk. We have not studied every 
possible situation, yet we have presented several models 
that enables us to describe more general situations. A real 
disk is thought to accrete matter both from its plane and 
outside it, with parameters varying with time. Thus, a real 
accretion is likely to be a succession of elementary events, 
each of them being represented by one of our models. To 
illustrate this, we now present the evolution of a disk that 
accretes gas in its plane over a Hubble time, and from out 
of its plane during 5 Gyrs at the middle of its life. 

In this run, the initial bulge to disk mass ratio is 0.7, 
and the halo to disk is 0.75. In-plane accretion is im- 
plemented as before with a rate that makes the initial 
disk mass double in 10 Gyrs. From t = 6 Gyr to t = 11 
Gyr, out-of-plane accretion takes place with R acc = 0.8, 
\P acc = 25, and a mass flow corresponding to the initial 
disk mass in 8 Gyrs. 

Before t = 6 Gyr, the bar evolution is rather simi- 
lar to the 2-D run A. Here, the initial bar is very weak- 
ened at t = 5 Gyr, while its strength had reached 0.68. 
Then, a new bar appears in the disk around t = 6 Gyr 
(see Fig. |l^). This bar is initially rather weak, then out- 
of-plane accretion strengthens it, as it has already been 
shown. The maximal bar strength is 0.92. In the out-of- 
plane runs with a maximal bar strength larger than 0.8, 
the bar strength drops below 0.4 in about 2 Gyrs. Yet, the 
very strong bar formed by out-of-plane accretion is here 
maintained by in-plane accretion. For instance, at t — 10.5 
Gyr, the bar has the same aspect as at t — 8 Gyr, as shown 
on Fig. [l9[ and is bar strength still is 0.85. Later, a bar 
whose strength is larger than 0.34 remains until the end 
of the simulation, at t = 15 Gyr. 

Concerning the evolution of the bar pattern speed 
(normalized to 1 at t = 0), the second bar formed around 
t = 6 Gyr is faster than the initial one (1.36), and the pat- 
tern speed is triggered by the out-of-plane event around 
t = 8 Gyr (2.35). After being very strong, the bar first 
slows down (1.97), but is then re-accelerated by the main- 
tenance effect, while it becomes shorter, as in the 2-D run 
B: the bar pattern speed at t = 15 Gyr is finally 2.31. 
Thus, the bar pattern speed doubles or even more over a 
Hubble time, as was the case for runs A and B (Sect. |3. 



5. Discussion of numerical models 



5.2 



5.3 



General discussions are given in subsections and 
A discussion of the ISM scheme is given in subsection 

5.1. 2-D models 



It is normal to expect that gas accreting galaxies main- 
tain a relatively high gas fraction in their disk. However, 
even if the disk mass doubles in a few Gyrs, the gas to 
total (stars +gas) mass ratio is generally not more than 
0.5: the star formation and gas consumption rate is of the 



same order as the accretion rate. The long-term evolution 
of a gas accreting disk depends not only on the dynam- 
ics of this gas, but also on the evolution of stars formed 
in accreted matter, particularly for the consequences on 
stellar bars. Mass exchanges between stellar and gaseous 
phases have to be treated realistically. Our simulations 
thus include detailed star formation and stellar mass-loss 
schemes. In particular, the stellar mass-loss is not approx- 
imated as instantaneous: this is rather important, since 
mass re-injection over several Gyrs would change radial 
flows of matter, which then influence the bar properties 
(see Jungwiert et al. 2001). However, there are still some 
shortcomings in the present model, since gas is assumed 
isothermal and the dependence of the yield on metallicity 
-that evolves with time, and also with accretion- is not in- 
cluded; it is hoped that these are only minor effects, that 
will not significantly affect the mass accretion rate. 

The dissipation time-scale k for the gas dynamics can- 
not be exactly justified, even if a realistic value is of a few 
Gyrs. Yet, this adjustable value has no major influence on 
the morphology of disks. When this time-scale k is var- 
ied between 5 and 10 Gyrs, the long-term evolution of the 
disk remains unchanged: the gas reaches an equilibrium 
between dissipation and heating by density waves which 
is not significantly influenced by the value of k. 

5.2. 3-D models 

Our three-dimensional simulations include all the major 
phenomena that influence the properties of bars: gas dy- 
namics, star formation, and stellar mass-loss and feedback. 
However, the resolution of our 3-D models is lower than 
in 2-D: the cell dimension is 450 pc (instead of 100 pc). 
The dark matter is represented by a rigid potential, and 
there are still many physical processes that are simplified 
or ignored (the multiphase nature of the gas, for instance). 
May be the most relevant uncertain part is the rate of star 
formation, that might not be a Schmidt law. 



5.3. Gas dynamics 

Different schemes were used to treat the gas dynamics 
in our simulations. They have given comparable results. 
They do not include a precise description of the real ISM. 
A more precise study may treat gathering and fragmen- 
tizing of gas clouds, and energy feedback from star forma- 
tion. However, a simplified description of the ISM seems 
to be suitable for a first approach of large-scale accretion. 
Actually, what may be more questionable is the dissipative 
scheme, for it does not conserve exactly angular momen- 
tum. Sticky-particles codes are generally thought to be 
rather convenient ; yet, we study gas- rich galaxies, which 
makes such a shortcoming be suspect. We have shown that 
angular momentum transfers and related mass exchanges 
are here fundamental processes, so we have to verify that 
the local non-conservation of the angular momentum does 
not invalidate our results. This is done in appendix |a[ 
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Fig. 19. Stellar and gaseous particles viewed perpendicular to the disk plane at t —6, 8, and 10.5 Gyr, for the in- and 
out-of-plane accretion run. 



6. Conclusion 

The numerical simulations presented in this work reveal 
how different can be the evolution of spiral galaxies with 
and without important gas accretion. While an isolated 
galaxy may spontaneously form a bar, consume its gas 
in stars, and when the stellar component is hot enough, 
see its bar weaken to a lens, and evolve definitely to an 
early-type system, a galaxy with significant gas accre- 
tion, can reverse this evolution, and experience several 
new bar episodes. The detailed processes are controlled 
by the strength of the bar and its gravity torques. The 
accreted gas may stay for a while in the outer parts of 
the disk, prevented to flow in by a strong bar; as soon as 
the bar weakens, the gas will progressively populate the 
disk, and form new stars there. This increases the self- 
gravity of the disk with respect to the stabilizing effect 
of the central bulge, favoring a new bar instability. Under 
the influence of this new bar, the material inside its coro- 
tation is driven inwards, and through vertical resonance, 
thickens to enrich the bulge again. This regulates the bar 
strength. 

In the course of the simulations, three to four bar 
episodes have been followed. An interesting discovery is 
that the pattern speed of the bar is increasing form one 
bar to the next. This can be considered as a rejuvenat- 
ing process. Indeed, the normal evolution of a bar is to 
slow down, while trapping more and more material, and 
lengthening, so that its corotation radius increases (a bar 
cannot extend beyond its corotation, since periodic orbits 
are perpendicular here). When another bar forms, it is 
shorter than the previous one. Its pattern speed is higher, 
not only because there is then more mass contained in a 
given radius, due to mass accretion. The pattern speed is 
even higher than what would be expected from the square 
root of mass increase. This increase may be related to 
the fact that mass is more and more concentrated during 
this evolution, and the galaxy is shifting progressively to 
early-types, with massive bulges. The bar is then shorter 
with respect to the galaxy radius (as shown by Combes & 
Elmcgrccn 1993). Also, gas accretion has the effect of in- 
creasing the inner angular momentum per unit mass. The 
bar is shorter and rounder. 

In the course of evolution with gas accretion, the ra- 
dius of the galaxy expands, a phenomenon also controlled 
by the bar strength, which gravity torques can maintain 
or expel the gas outwards. One of the most important con- 
sequence of gas accretion, is to maintain an almost perma- 
nent (though recurrent) spiral structure through the disk, 
and also a significant star formation rate along a Hubble 
time. This solves the problem of explaining the almost 
constant star formation rate of spiral galaxies along their 
history (e.g. Kennicutt 1983), and also the presence until 
today of conspicuous spiral structure in most spiral disks. 



The possibility to accrete gas from a direction highly 
inclined with respect to the galaxy plane introduces new 
processes that can be also favorable to bar reformation. 
The gas may not have to wait for a weakening of the pre- 
vious bar to flow in, and the distribution is then asymmet- 
ric, with an m=2 extra contribution, which makes bars 
become much stronger. 

One of the consequences of the evolution described in 
this work is that the bar pattern speed can no longer give 
the evolution state or " age" of a barred galaxy. Significant 
gas accretion can reverse the evolution and increase the 
pattern speed. Also, the argument that dynamical fric- 
tion would have slowed down the bars along their evolu- 
tion would have to be revised (e.g. Debattista & Sellwood 
1998). Another consequence is that the morphology of a 
galaxy, according to its environment, and the availability 
to accrete cold gas, may oscillate from late to early- types 
three or four times, before inexorably shift to an early-type 
galaxy. 

The maintenance of conspicuous spiral structure in 
most disk galaxies may mean that galaxies have sufficient 
gas to accrete in their outer parts, to constantly renew 
disk instability, even in the presence of a heating stellar 
disk and increasing bulge mass. The relative absence of 
bars in galaxies at high redshift (e.g. van den Bergh et al. 
1996) could be interpreted as a very short life-time of the 
bar episode, in the presence of dominating gas fraction in 
very young galaxies. Indeed, the fraction of barred galax- 
ies at a given epoch is interpreted, in this recurrent bar 
scenario, as the relative time spent by a galaxy in the bar 
episode relative to the un-barred one. 

The present work has only considered general prop- 
erties of gas accreting spiral galaxies, but many parame- 
ters (mass distribution, morphological types, geometry of 
the accretion, etc..) have to be explored before all conse- 
quences of accretion on galaxies can be drawn, which will 
be done in future work. 

Appendix A: Numerical schemes for gas dynamics 
and non conservation of integrals of 
motion 

A .1. Variations of angular momentum and energy 

We have shown that redistribution of angular momentum 
inside the disk by spirals and bars waves was partly re- 
sponsible for an increase in the bar pattern speed, which 
is an important consequence of gas accretion. However, 
the numerical codes we have used do not conserve an- 
gular momentum. Both ISM models, viscosity scheme or 
sticky-particles code, dissipate the angular momentum of 
gas clouds, at least with the chosen parameters. This could 
be a real process if a spin was stored in gas particles, but it 
is not the case in our simulations. Thus, we discuss in this 
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Integral of motion 


Value after 5 Gyrs 


2-D 


viscosity scheme 






Angular momentum 


0.97 




Kinetic energy 


0.60 


3-D 


sticky-particles scheme 






Angular momentum 


0.98 




Kinetic energy 


0.53 



Table A.l. Relative variations of angular momentum and 
energy of gas particles over 5 Gyrs, for the control run 
without accretion, star formation, and stellar mass-loss. 
The initial values are normalized to 1. 



appendix the non-conservations of angular momentum by 
gas dynamics schemes, in relation with our results, and 
more generally in all N-body simulations. Another basic 
integral, the energy, is not conserved in the ISM: here it 
corresponds to a physical process, for gas clouds radiate 
away a significant energy. 

In a sticky-particles code, when two gas particles col- 
lide, their local angular momentum is generally not con- 
served. According to equations |is| to 24, the angular mo- 
mentum of the particles in the referential of their mass 
center before the collision is: 



(r 2 - ri) x v t 



mim,\ 
M 



(A.l) 



After the collision, only velocities have changed, and the 
angular momentum is: 



L' = (r 2 - ri) x v' t 
L' = (3 t L 



m2TOi 
M 



(A.2) 
(A.3) 



Thus, the sticky-particles code conserves the angular mo- 
mentum only when /3 t = 1. It does not when /3 t = 
0.65. However, during a collision, the angular momentum 
around the mass center is changed, but the momentum 
of the mass center around the galactic center remains un- 
changed. Actually, the second one is much larger: its mean 
value, for each pair of gas particles susceptible to collide, is 
almost 10 3 larger than the mean value of the first one (this 
measure is for run 0F1). Thus, only a small fraction of the 
angular momentum can be reduced. We are now to show 
quantitatively that the variations of angular momentum 
induced by both dissipative schemes are not as important 
as the effects of gas accretion that we have shown. 

Accretion is responsible for changes in angular momen- 
tum and energy, thus, in order to quantify the non conser- 
vation of integral of motions due to our code, we compute 
the variation of energy and angular momentum in an iso- 
lated disk (run 0F1 in Sect. ||). Star formation and stellar 
mass-loss modify the gas content and the importance of 



gas dynamics, thus we also suppress them. In table A.l 
we give the relative variation of angular momentum and 
energy in the gaseous disk over 5 Gyrs. The quantities 
taken away by particles going out of the grid have been 
accounted for, thus the gas dissipation is the only cause for 
the observed variations. The ISM undergoes an important 




3 4 5 
time (Gyr) 

Fig. A.l. Bar strength in a 2-D simulation. Solid line: 
angular momentum is conserved (sticky-particles (3t = 1), 
dashed line: angular momentum is not conserved (viscos- 
ity scheme), dotted line: angular momentum is conserved 
(modified viscosity scheme). We show the scatter (error 
bar) measured on 5 simulations starting with 5 different, 
randomly chosen, initial position of the particles, all the 
physical parameters remaining unchanged. 

loss of energy, that is physically justified, and a small loss 
of angular momentum: as explained qualitatively before, 
collisions in the 3-D code, and viscosity in the 2-D one, 
dissipate only a small fraction of the angular momentum, 
that cannot be compared with the large supply of angular 
momentum by accretion. 

A.2. Influence of the angular momentum conservation 
on previous results 

We have explained that numerical schemes for gas dy- 
namics do not conserve exactly the angular momentum. 
We wish to show that this is not a critical point, and does 
not weak our conclusions about gas accretion. 

The dissipative schemes we have used suppress local 
vortices, and tend to reduce the angular momentum of 
two colliding particles. This may modify the dynamics of 
the galactic center. Thus, we modify the value of (3 t to 1, 
then the angular momentum is conserved. We test sim- 
ulations with this parameter against previous results es- 
tablished with a sticky-particles code, i.e. 3-D results. For 
2-D results, we compare the viscosity scheme to a modified 
version in which only the radial velocities are reduced, still 
in order to conserve angular momentum: in this modified 
version, a particle at a position r from the disk center, 
with a velocity «, in a cell in which the mean velocity is 
< v >, will undergo a force that reads: 



/ = -k {{v- < v >) • r) ■ 



(A.4) 



We also compare the 2-D viscosity scheme with a 2-D 
sticky-partic les c ode, using t = 1. 

On Fig. AA, we compare the evolution of the bar 
strength in a 2-D disk for the viscosity scheme, its mod- 
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Fig. A. 2. Bar pattern speed in a 3-D simulation. Solid 
line: angular momentum is conserved (flt = 1), dashed 
line: angular momentum is not conserved (j3t = 0.65). The 



error bar has the same meaning as on figure AT 



ified version that conserves angular momentum, and the 
sticky-particles code with (3 t = 1- Fig. [A.2| shows the evo- 
lution of the bar pattern speed in a 3-D run, for the sticky- 
particles code with j3t = 0.65 and 1. The results are quite 
similar for the different codes, and the differences between 
each of them is not significant: they do not largely exceed 
the scatters obtained on 5 simulations starting with the 
same values for physical parameters, but with 5 different, 
randomly chosen, initial positions for the particles. 

For both 2-D and 3-D models, we conclude that con- 
serving the local angular momentum or not does not sig- 
nificantly influence the bar evolution and the spiral struc- 
ture. This confirms all our previous results. More gener- 
ally, this shows that the 2-D linear dissipative scheme and 
the sticky-particles code are rather convenient when one 
studies the evolution of disk galaxies. Furthermore, the 
dissipation in the real ISM is likely to be somewhere in be- 
tween our two schemes. One scheme exactly conserves the 
angular momentum, the other one annihilates local vor- 
tices: real gas clouds may dissipate a fraction of angular 
momentum in cloud-cloud collisions, and store it in their 
spin. As the results are similar with or without dissipat- 
ing angular momentum, it seems likely that the response 
of the real ISM to gas accretion on galactic disk is similar 
to what we have drawn out from our simulations. 
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